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Meson decays in which neither constituent quark ap- 
pears in the final state are expected to be dominated 
by annihilation diagrams, in which the two quarks in- 
teract directly. Such processes provide interesting in- 
sights into the internal dynamics of B mesons and need 
to be understood to make precise predictions on B me- 
son decays. Such diagrams cannot be calculated by as- 
suming factorization since both the quarks play a role 
and a reliable theoretical prediction for the correspond- 
ing amplitudes does not exist. These amplitudes are 
expected to be suppressed with respect to amplitudes 
where one of the two quarks is a spectator by a factor 
- f s /m B - 0.04 (f fl - 200 MeV and m B = 5.28GeV/c 2 
are the B meson decay constant Q and mass, respec- 
tively). This factor represents the amplitude for the two 
quark wave functions overlapping, a necessary condition 
in annihilations. So far no process relying entirely on 
annihilation has been observed and the assumption that 
these types of diagrams can be neglected is frequently 
used in theoretical calculations. Some studies [2( indi- 
cate, though, that processes with a spectator quark can 
contribute to annihilation-mediated decays by rescatter- 
ing if the final state is reached in two steps: a decay into 
two mesons that can occur with a spectator quark, and 
a subsequent strong interaction between the two mesons 
which produces the final state of interest. Figurc^shows 
the Feynman diagram for the decays B + — yD^ + K^ and 
B + ^Df-ir° Q, and the hadron-level diagram for the 
rescattering. Strong rescattering could then mimic large 
contributions from annihilation diagrams to the level of 
not being negligible any more. 

The decays B + -^D^ + K^ are particularly suited to 
study annihilations because of their relatively clean ex- 
perimental signature and because their branching frac- 
tions are expected to be at the level of the current sensi- 
tivity (10 -5 ) if large rescattering occurs, or three orders 
of magnitude below if not Q. Moreover the branching 
fraction of these decays can be used to constrain the anni- 
hilation amplitudes in the phenomenological fits Q that 
allows to translate the measurement of the amplitude 
of B + ^D°K + into estimates of the Cabibbo-suppressed 
decay B°^D°K° needed in some CP measurements 
Neither of the modes studied here has been observed so 
far, and a 90% confidence level upper limit on the branch- 
ing fraction B{B+^D*+K°) < 9.5 x 10~ 5 has been es- 
tablished by CLEO @. 



In this paper we present the results of the search for 
B+_ -► D^+K° decays in 225.9 ± 2.5 million JUS) -> 
BB decays, collected with the BABAR detector at the 
PEP-II asymmetric-energy B factory at SLAC. We use 
a Monte Carlo (MC) simulation of the BABAR detec- 
tor based on GEANT4 Q to validate the analysis pro- 
cedure, estimate efficiencies, and to study the relevant 
backgrounds. We also use 12.4 fb of data collected at 
a center-of-mass energy approximately 40 MeV below the 
T(45) mass. 

Candidates for D + mesons are reconstructed in the 
modes D + — > K~it + ir + and D + — > K®ir + . Candidates 
for D* + mesons are reconstructed in the mode D* + — > 
D°tt + , where the D° subsequently decays to one of the 
four modes K~ir + , K~ir + Tr Q , K~ir + ir~ir + , or K®tt +/ k~. 

K® candidates are reconstructed from two oppositely- 
charged tracks with an invariant mass 491 < m^+ n - < 
504 MeV/ c 2 (corresponding to a ±2 standard deviations, 
a, window around the mean value in control samples). 
The x 2 °f the 7r + 7r~ vertex fit must have a probability 
greater than 0.1% and the K® flight distance from the 
primary vertex in the plane transverse to the beam axis 
in the event must be greater than 2 mm. Kaons and pi- 
ons coming from the D are required to have momentum 
in the laboratory frame greater than 200 MeV/c and 150 
MeV/c, respectively, except in the decays D° — > K~tt + 
(K~ir + ir°) where the momentum threshold for both the 
tracks is 200 (150) MeV/c. To identify charged kaons 
we use a selection with an efficiency of 95% and a 12% 
pion misidentification probability. 7r° candidates are re- 
constructed combining two photons with invariant mass 
120 < TO 77 < 150 MeV/c 2 (corresponding to a ±2<r win- 
dow around the mean value estimated on control sam- 
ples) and a minimum total energy in the laboratory frame 
of 200 MeV. For the D° -> R-tt+tt decay we select the 
dominant resonant contributions with a requirement on 
the Dalitz density distribution [9j. 

Finally, the D + and D° candidates are required to have 
an invariant mass within 2a of the mean values. The D + 
and D mass resolutions are mode- dependent and range 
between 5 and 8 MeV/c 2 . We form D* + candidates by 
combining D° candidates with charged tracks. The mass 
difference between the D* + and the D° candidates is re- 
quired to be within 2a of the mean value as estimated 
on control samples. The resolution is mode-dependent, 
approximately 0.6 MeV/ c 2 in all cases. We combine D + 
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FIG. 1: Annihilation diagram for the decay B+~^D M+ K° s (\ett), tree diagram for B + — *Ds 7r°(center), and hadron-level 
diagram for a possible rescattering contribution to B + —*D^ + Kg via 



or D* + candidates with a K® to form B + candidates. To 
improve the resolution on the four-momentum of all the 
intermediate composite particles we apply a kinematic 
fitting technique that constrains their masses to the nom- 
inal value |l(| and their charged daughters to come from 
the same vertex. 

We only accept events with a reconstructed candidate 
and a total measured energy greater than 4.5 GeV, de- 
termined using all charged tracks and neutral clusters 
in the electromagnetic calorimeter with energy above 30 
MeV. The remaining background comes predominantly 
from continuum qq production. This background is sup- 
pressed using variables that characterize the topology of 
the event. We require the ratio of the second and ze- 
roth order Fox- Wolfram moments ^lj to be less than 
0.5. We compute the angle 6t between the thrust axis 
of the B-meson candidate and the thrust axis of the rest 
of the event. The thrust axis is defined as the direction 
that maximizes the sum of the longitudinal momenta of 
the particles in the center-of-mass (cm.) frame. In this 
frame BB pairs are produced approximately at rest and 
have a uniform |cos#t| distribution. In contrast, qq pairs 
are produced back- to-back, which results in a |cos#t| dis- 
tribution that peaks at unity. To further suppress back- 
grounds we use a Fisher discriminant T constructed from 
the scalar sum of the cm. momenta of all tracks and pho- 
tons, excluding the B candidate decay products, flowing 
into nine concentric cones centered on the thrust axis 
of the B candidate ^j. The more spherical the event, 
the lower the value of T. Figure |21 shows the distribu- 
tion of J- and |cos#t| on signal MC and on off-resonance 
data, which contain exclusively continuum qq production 
events. We also exploit the charge correlation between 
the B and the leptons and kaons produced in its decays 
to classify the events in three mutually exclusive cate- 
gories with different levels of contamination from contin- 
uum background: events with at least one lepton with 
charge opposite to the B candidate, events with no lep- 
ton and at least one kaon among the tracks that do not 
form the B candidate but have opposite charge, and all 
the other events. The optimization of the selection is 
performed separately for each decay mode and for the 
three categories by maximizing the ratio of signal effi- 



ciency, estimated with MC, over the square-root of the 
expected number of background events, estimated in data 
sidebands: the maximum allowed value of |cos#t| ranges 
between 0.8 to 1 (i.e. no cut) and the maximum allowed 
value for T varies from 0.1 to 0.7. 



0.05 - 





1 

{[ 






{ 






n 

\ 


1 

I 


1 


T 


U - 




T 






• 


vv. 


, ,-rf <m 1 





2 

Fisher(a.u.) 




0.25 0.5 0.75 1 
|cos6 T | 



FIG. 2: Distribution of the discriminating variables |cos#t| 
and T in the B + K° signal MC (histograms) and the 
off-resonance data (dots). 



mES 



We extract the signal using the kinematic variables 



■>*2 



, where E^ is the beam energy in the cm. frame, p* is 
the cm. momentum of daughter particle i of the B meson 
candidate, and rrii is the mass hypothesis for particle i. 
For signal events, ttt-es peaks at the B meson mass with a 
resolution of about 2.5 MeV/c 2 and AE peaks near zero, 
indicating that the B candidate's total energy is consis- 
tent with the beam energy in the cm. frame. The AE 
signal band is defined as \ AE\ < 2.5a and within it we de- 
fine the signal region as 5.2725 < tues < 5.2875 GeV/c 2 
and the ibes sideband region as 5.2000 < ities < 5.2725 
GeV/c 2 . The AE resolution a is mode- dependent and 
approximately 18 MeV. We also define the AE sideband 
region as 2.5cr < \AE\ < 0.12 GeV and 5.2 < m E s < 5.3 
GeV/c 2 . Table |U shows the efficiency for each sub-decay 
mode estimated with simulated events. Depending on 
the mode, in 1.5 to 7% of the events there is more than 
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TABLE I: Efficiencies for the B + -» D w+ K° candidate re- 
construction in each sub-decay mode. The branching fraction 
of the Z)(*) + decay chains considered [Til ] are also shown. 



D mode 


£i(%) 


B(%) 


D+ -f 


K U S TT + ; Kg -» 7T-7T + 




17.3 


0.97 ± 0.06 


D+ 


if~7T + 7r + 




17.7 


9.2 ± 0.6 










D° 


-> if~7T + 




18.5 


2.57 ± 0.06 


D° 






6.4 


8.8 ± 0.5 


£>° 


— ► if ~7r + 7T _ 7r + 




10.1 


5.05 ± 0.21 


D° 




"7T+ 


10.6 


1.37 ± 0.08 




one B candidate. We select the B candidate whose 
candidate's mass is closest to its nominal mass or, in case 
two B candidates are formed by the same the flW can- 
didate, one with the smallest value of |A23|. 

After the selection described above, two classes of 
backgrounds remain. First, there is combinatorial back- 
ground in the signal region, coming from random com- 
binations of tracks in the event. We estimate this 
background from the sideband of the tues distribution, 
describing it with a threshold function dN/dmES °c 

toes\A - m Es/ E b 2 ex P [ - £ (! ~ m !ks/ E b 2 )] > charac- 
terized by the shape parameter £ [13J. We obtain the 
parameter £ from a fit to the distributions of tues m 
data, in the AE sideband region. The number of com- 
binatorial background events is obtained by scaling the 
events in the sideband of the tues distribution into the 
signal region with the ratio of the threshold function area 
in the two regions. Including systematic errors, we esti- 
mate 56.3 ± 3.0 and 22.0 ± 1.8 events for the B+^D + K° S 
and B + ^>D* + K® mode, respectively. Second, there is 
peaking background due to misreconstructed B meson 
decays that have an m,ES distribution peaking near the 
B mass. We study the peaking background with MC and 
we estimate it to be 4.4 ± 1.2 and 1.2 ± 0.6 events for the 
B+^D+K° S and B+^D*+K° modes, respectively. The 
dominant contribution to the peaking background comes 
from well-known B° -> D^~X + decays (X+ = tt+, p+, 
af). As a cross-check, we also estimate the peaking back- 
ground using candidates from the D mass sidebands in 
data and we find results consistent with the MC predic- 
tion. 

Figure|21shows the m,Es distributions in the AE signal 
band for the two modes after the selection. The expected 
background is superimposed. 

To compute the confidence level (C.L.) at which the 
data agree with a given hypothesis on B(B+^D^ + K°) 
we use a frequentist technique [l4|, which treats prop- 
erly the small number of events and includes the sys- 
tematic errors directly in the computation of confidence 
intervals or limits. The C.L. is defined as the fraction 
of times a random number, following the expected dis- 



FIG. 3: The m E s distribution for the a) B + -*D + K° S and b) 
B + ^D* + K® candidates within the AE signal band in data 
after all selection requirements. Combinatorial (full line) and 
peaking (dashed line) backgrounds are superimposed. 



tribution of the number of events in the signal region 
(N exp ), exceeds the number of observed events (N canc t 
in Tab. ITT|l . N exp is distributed according to the sum of 
Poissonian distributions with mean values /i distributed 
as follows: for a given value of B{B+^D^+K°) we esti- 
mate [i as the sum of the expectation value of the number 
of events from the combinatorial and peaking background 
( N como and N pea k, respectively), and from the signal 

{N s ig), fl — N cam b + Np ea k + N sig . 

We estimate N com b by scaling the number of events 
in the niES sideband to the signal region and by con- 
sidering the Poisson fluctuations of the number of events 
in the sideband and the systematic uncertainties on the 
threshold parameter £. We estimate N pea k from the MC, 
taking into account its limited statistics. TablelTTlreDorts 
the mean values and standard deviations for N comb and 
Np e ak- Finally, for a given value of the branching frac- 
tion, N S ig is obtained as: 

N sig = B(B+ -> D&+K ) xN b x EidBi (1) 

where the number of B^ mesons (Ng) and the prod- 
uct of the efficiency and the branching fraction of the 

TABLE II: The number of candidates in the signal region in 
data (-/V can d), the corresponding expected combinatorial back- 
ground (iV com b), the peaking background (iVpeak), the prob- 
ability (Pbkgd) of the data being consistent with the back- 
ground fluctuating up to the level of the data in absence of 
signal, and the 90% confidence- level upper limit. Systematic 
uncertainties are included. 



B mode 


A'cand 


■^Vi comb 


A'poak 


Pbkgd(%) 


90% C.L. 


D + K» 


57 


56.3 ± 3.0 


4.4 ± 1.2 


69 


0.5 x 10" 5 


D*+ K» 


28 


22.0 ± 1.8 


1.2 ± 0.6 


24 


0.9 x 10~ 5 
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sub-decay modes (ZJiCiBi) are varied according to Gaus- 
sian distributions within their systematic uncertainties. 
The systematic errors on the reconstruction efficiency are 
shown in Table HTTl and include the uncertainty due to lim- 
ited MC statistics, uncertainty on tracking efficiency, K° 
and 7T° reconstruction, charged-kaon identification, other 
selection criteria. They have all been estimated by com- 
paring the data and simulation performances in control 
samples. Also, the uncertainties on Nb (1.1%) and on 
the branching fraction of the sub-decay modes have been 
taken into account. The total uncertainty is obtained by 
adding the contributions from the individual sources in 
quadrature. 

Calculating the C.L. with the procedure just described 
and setting B{B+—>D ( -*'> + K°) = 0, we estimate the 
probability of the background to fluctuate above the 
observed number of events to be 69% and 24% for the 
B+^D+K° and for the B+^D*+K° modes, respec- 
tively. In absence of significant signal we then set the 
following upper limits on the values of the branching 
fractions corresponding to a C.L. of 90%: 

B(B+^D+K°) < 0.5 x 10^ 5 , (2) 
B(B + ^D*+K°) < 0.9 x 10" 5 . 

We also compute the branching fractions 
B(B+^D+K°) = (-0.28lg;|J) x lO" 5 and 
B(B+^D*+K Q ) = (0.28±{J3J) x lO" 5 . The errors 
above include both the statistical and the systematic 
uncertainties. 

In conclusion, we report on the search for the rare 
decays B + — >D(*' + K®, which are predicted to proceed 
through annihilation diagrams. We do not observe any 
significant signal and we set 90% C.L. upper limits on 
their branching fractions. 
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TABLE III: Relative systematic errors on the branching fraction due to, respectively: MC statistics , track reconstruction , 
Kaon identification, K% and n reconstruction efficiencies, and the data-MC agreement on the signal shapes of AE, cosd T , and 
T. 



D mode 


MC(%) 


Tracks(%) 


Kaon(%) 






AE(%) 


cos(9 T (%) 




Total(%) 


D+ -> K»n+ 


1.4 


1.1 




0.3 




0.4 


0.4 


0.7 


2.0 


D+ -f K--K+-K+ 


0.5 


1.1 


0.4 


0.2 




0.4 


0.4 


0.7 


1.5 


- K-7V+ 


0.9 


1.2 


0.4 


0.2 




0.8 


0.4 


0.7 


1.8 




0.3 


0.4 


0.1 


0.1 


0.3 


0.2 


0.1 


0.3 


0.7 




0.5 


0.9 


0.2 


0.1 




0.1 


0.2 


0.4 


1.2 


D° - K°7r+7r- 


1.0 


1.0 




0.2 







0.2 


0.4 


1.5 



